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The advancement of fabrication and lithography techniques of semicon- 
ductors have made it possible to study bi-layer systems made of two electronic 
layers separated by distances of several hundred Angstroms. In this situation 



'^ ' the electrons in layer 1 are distinguishable from those in layer 2, and can com- 

Q ' municate through the direct inter-layer Coulomb interaction. In particular, 

if a current is applied to one of the layers, the electrons in the second will 

>' 
("j^ ' be dragged giving rise to a transresistance pD- In this article we review re- 

^n . 

^^ ' cent theoretical and experimental developments in the understanding of this 

Qv ' effect. At very low temperatures it turns out that phonons dominate the 

transresistance. The direct Coulomb interaction and plasmon excitations are 

2 . 

important at temperatures T > O.lTp, with Tp the Fermi temperature. If a 

^ ' magnetic field is applied the transresistance is increased, in a very interesting 

o' 

O ■ interplay between po and Landau quantization. The non-dissipative drag is 



also reviewed. 



I. INTRODUCTION 



Electron-electron (e-e) interactions are responsible for a multitude of fascinating effects 
in condensed matter. They play a leading role in phenomena ranging from high temperature 
superconductivity and the fractional quantum Hall effect, to Wigner crystalization, the 
Mott transition and Coulomb gaps in disordered systems. The effects of this interaction 
on transport properties, however, are difficult to measure. A new technique has recently 
proven effective in measuring the scattering rates due to the Coulomb interaction directlya. 

This technique is based on an earlier proposal by PogrebinskiiQ'B. The prediction was 
that for two conducting systems separated by an insulator (a semiconductor-insulator- 
semiconductor layer structure in particular) there will be a drag of carriers in one film due 
to the direct Coulomb interaction with the carriers in the other film. If layer 2 is an "open 
circuit", and a current starts fiowing in layer 1, there will be momentum transfer to layer 2 
that will start sweeping carriers to one end of the sample, and inducing a charge imbalance 
across the film. The charge will continue to accumulate until the force of the resulting 
electric field balances the frictional force of the interlayer scattering. In the stationary state 
there will be an induced, or drag voltage Vo in layer 2. 

A fundamental difference emerges in the configuration of the current drag experiment 
with that in which the in-plane resistance is measured. For a perfectly pure, translationally 
invariant system, the Coulomb interaction cannot give rise to resistance since the total 
current commutes with the Hamiltonian H. This means that states with a finite current 
are stationary states of H and will never decay, since the e-e interaction conserves not only 
the total momentum but also the total current. (For electrons moving in a periodic lattice, 
momentum and velocity are no longer proportional and the current could in principle decay 
by the e-e interaction.) If the layers are coupled by the Coulomb interaction, the stationary 
states correspond to linear superposition of states in which the current is shared in different 



amounts between layers: the total current within a given layer is not conserved and can 
relax via the inter-layer interaction. 

This mechanism of current degrading was studied in the pioneering experiment of 
Gramila et a/.a for GaAs layers embedded in AlGaAs heterostructures. The separation 
between the layers was in the range 200-500A. The coupling of electrons and holesiil and the 
coupling between a two dimensional and a three dimensional was also examined^. 

If we call / the current circulating in layer 1, the drag resistance (or transresistance) is 
defined as 

The naive expectation of the temperature dependece of p^ is that it should vanish as T^ 
at low temperatures. This results from the exclusion principle which limits the scattering 
to states within ksT of the Fermi surface, and there is one factor of T per layer. The first 
experiment seemed compatible with the e-e mechanism with some discrepancy since the 
observed ratio Po/T"^ was not a constant at low temperatures. This experiment motivated 
a rather extensive theoretical effort to understand in better detail the mechanisms of cross- 
talking between layers between which electrons are not allowed to tunnel. The theory 
in turn stimulated new experiments and in the recent years we have seen Coulomb Drag 
developing as a field in itself. In the present article we briefly review the progress made 
in understanding different aspects of the problem. Section |I| provides some background on 



the Coulomb interaction and screening for double layer systems. In Section |T| we discuss 
the perturbative equation for pj^ used in most treatments. This equation can be derived in 
various ways. In the Appendix we present the memory function derivation and the detailed 
low temperature dependence assuming Coulomb scattering only. 

The following sections discuss the crossover between three regimes. Phonon exchange, 
discussed in Section 0, dominates the drag at T < O.lTp, with Tp the Fermi temperature. 
For T K, O.lTp , single particle Coulomb scattering is the dominant effect, and for T > 0.2Tp, 
plasmons (See section [V|) are responsible for an enhancement of the drag current. The 



theoretical studies of the effects of disorder and locahzation are reviewed in Section 0. The 
enhancement of pd due to an apphed magnetic field perpendicular to the layers and the 
interplay between Landau quantization and interlayer scattering is reviewed in Section |Vflj . 
Finally in Section |V111| we discuss the theory of non-dissipative drag and some experimental 



attempts to measure it. 

II. INTER-LAYER COULOMB INTERACTION 

The leading actor in this play is the Coulomb interaction between electrons in different 
layers. This interaction is responsible for the scattering, and will be screened by the density 
fluctuations within each layer. If one neglects the interaction between layers, the effective 
interaction V{q,uj) within a given layer calculated in the Random Phase Approximation 
(RPA) is given byQ 

1 + H(g)xw,^) 
with Vh{q) = 2Txe^ /q the bare interaction and x{Qi ^) the function characterizing the response 
of the charge 6p{q,uj) to an external potential ip{q,u!): 6p{q,uj) = —x{q,uj)ip{q,uj). For two 
coupled layers we also have the bare interlayer interaction Ub{q) = 27r e^ exp{—qd) /q, with 
d the distance between layers. In addition, the response functions could be different if the 
layers are not identical. The effective interaction can be written now as a 2 x 2 matrixpBcl 

V{q,u;) = Vs{Q) ^ (2) 

l + XW,^)Vij(g) 

with [VB(g)]ii = [VB{q)h2 = Vb{q), [VB{q)]i2 = [V^B(g)]2i = Ub{q), and [x{q,^)]i, = 
SijXiiq,^), with i = 1,2 labeling the layer. Note that xi and X2 could be different, so 
Eq. @ is valid for non-identical layers. The off-diagonal element of V{q,uj) is the screened 
interlayer interaction 

^^^^'''^^'' ^ [1 + Xi(g, oo)V,{q)] [1 + X2{q, uj)V,{q)] - U^,{q)Xi{q, ^)X2(g, ^) ^ ^^''^^^''^^ 

(3) 



The zeroes of the denominator in the above equation correspond to the collective modes 
of charge oscillations. It is illustrative to write the screened interaction for identical layers 
in the form (omitting the argument q in the interactions for simplicity) 



U''\q.^) = \ 



(4) 



1 + x(g, uj)iy, + Ub) 1 + x(g,o;)(H - U, 

which puts in evidence the symmetric and antisymmetric modes. The response function 
x(g, uo) has real and imaginany parts. The imaginary part corresponds to the continuum 
of particle-hole excitations that exist for frequencies uo < vpq- The collective modes corre- 
sponding to the poles of U^^^q^oj) will be non-decaying at T = if they fall out of the 
particle-hole continuum, which guarantees that the poles are in the real axis. We can find 
these stable modes by expanding the real part of x(?) ^) fo^' the free electron gas for small 
wave vectors and u ^ vpq: 

x(g,^)^ --(-)' (5) 

with m the electron mass, and n the electron densitytO. In the small q limit, V^+f/fe ~ Aiie^ /q^ 
and H — f/f, ~ 2TTe'^d. If we substitute this in Eq. (^ we find that the two poles are at 
frequencies u}±{q) given by 



V m 



27rnd 
V m 

This in-phase and out-of-phase modes are usually labeled "optical plasmon" and "acous- 
tic plasmon" respectively, and are stable at T = 0. At finite temperatures x(?) ^) acquires 
a finite imaginary part outside of the particle-hole continuum, the poles move into the com- 



plex plane and the plasmons acquire a finite lifetime. In section ^ we will discuss the role of 



plasmons in enhancing the drag resistance at temperatures of the order of the Fermi energy. 



III. FRICTIONAL DRAG 

The cross resistance is calculated in most of the recent theoretical treatments from the 
following expression, valid when the interlayer interaction is treated perturbatelyliiJ(See Fig- 
ure 1): 

In Appendix A we reproduce the derivation by of this equation by Zheng and MacDonald 
who used the memory function formalismllHl. This derivation is general and contains Equation 
(P) as the limit of lowest order in the inter-layer interaction. We will show here a more 
restricted derivation using the scattering formalism, which has the advantage of providing 
us with an interpretation of the origin of each of the terms in Equation (|^). The derivation is 
very close to the one presented in the paper by Jauho and SmithEla. Alternative treatments 
are the collective excitation approaches, and the Green's function formalism staring from 
the Kubo formulall^Ea. In the linearized Boltzman (or scattering) treatment, the spectral 
functions Xj(g,a;) are the susceptibilities of the non-interacting case. 



We compute first the rate of change of the momentum of the electrons in layer 2 due to 
the scattering from electrons in layer 1: 

- /ki+q(l - /ki)/k2-q(l - /ka)} ^(eki + ek2 " Ck^+q " Ck^-q)- (7) 

In the above equation, /k refers to the distribution of electron states in layer 1, and /^ 
to that of the electrons in layer 2. Since there is no current in layer 2, the corresponding 
electron distribution is the unperturbed, free-fermion equilibrium distribution. The electron 
distribution on layer 1, on the other hand, corresponds to a Fermi distribution displaced in 
fc-space by an amount mvi/h, with Vi being the drift velocity of the electrons in layer 1. 



The above expression can be interpreted as a sum of procesess according to Fermi's golden 
rule: the factor |f/(g)p is the square of the matrix element of the transition that involves a 
momentum transfer hc[, the various factors of / correspond to the probabilities of transitions 
from occupied to empty states, and the delta function ensures conservation of energy in the 
scattering process. Note that if layer 1 were in equilibrium, the magnitude in curly brackets 
is zero due to the detailed balance condition 

/k°(l - /k°')/k°"(l - /k°'") = /kVl - /u°)/k°'"(l - /k°"), (8) 

which is an identity if ek + ek' + ek" + ek'" = 0. There is therefore a finite momentum transfer 
to layer 2 due to the asymmetry of the electron distribution in layer 1. 
We now linearize in the electron distribution of layer 1 : 

/k = fLm../n = f^ - If ^k ■ VI = /° - ^/°(1 - f^)hk . VI. (9) 

Subtituting the linearized expressions for /i in Equation (^, and using the detailed 
balance condition the rate of momentum transfer has the form 

f - i^ / (If<''"'^|''<''iV (If / (Sp^td - /.%,)/£,(! - fi-.)H^. + ^., - ^..« - ^..-, 

(10) 

A couple of simple technical manipulations are needed to recast this expression in the 
form of our basic equation. First we use 

5(eki + Cka - eki+q - ek2-q) = / Muj5{huj - Ck^ + eki+q)5(^u; + e^^ - eka-q), (11) 

and 

/°(ek) [l - /°(ek + hu)] = [/0(ek) - f{e^ + hu;)] / [1 - exp{-hu;/kBT)] . (12) 

We then introduce the function x^il,^)'- 

(13) 



and substitute in (|T0]): 



dP2_v, 1 f ^q^(;,,)^|i.(,)p/rf^I^Xifa^)Imx^(g,c.)^ ^^^^ 



rft 2 ksT J (27r)2' ^' ' '^" J smh\hu/kBT) 

where we have omitted the vectorial nature of P2 and V2 since they are in the same direction. 
In order to compute the resistance, we need to equate the rate of momentum transfer to the 
total force per particle on the electrons in layer 2 due to the generated electric field E2: 

^ = +n2eE2, (15) 

and we get the final expression using p^ = E2/J1 with ji = riievi. 

The minus sign in Eq. (P) means that the induced drag voltage is opposite to the 
resistive voltage drop in the current carrying layer. This is so because the Coulomb-induced 
scattering sweeps the carriers along the dragged layer in the same direction as those in 
the drive layer. From the above equation we see that the drag resistance is a convolution 
of the density fluctuations within each layer, which at low temperatures are restricted to 
low frequencies by the factor sinh {Phuj/2) in the denominator. From the structure of 
the above equation we see that measurements of pD can provide information on the inter- 
layer scattering mechanism as well as the in-plane fluctuations, the information on which is 
included in x(g, cj). The low temperature prediction of Coulomb scattering alone isEJ 

_ m nCmkBTf 

^'' nenQhEFiqTFdyikFdy ^^ 

where gxF is the single-layer Thomas-Fermi screening wavevector, Ep is the Fermi energy, d 
is the interlayer spacing and kp is the Fermi wave vector. Due to the finite separation of the 
layers, the scattering is limited to small anglescJ. Large angle scattering events require large 
values of the momentum transfer q, and these processes are suppressed by the exponential 
dependence ~ e~'^'^ in the Fourier transform of the inter-layer Coulomb interaction. In 
two dimensions, the phase space for backscattering diverges and gives rise to logarithmic 
corrections to the simple T^ dependencec3. This corrections are missing here due to the 
suppresion of backscattering. The screening is also enhanced at small wave vectors, and is 



more effective as the separation between layers increases [see Eq (pll)]. The combination of 
the suppression of backscattering and the enhanced screening gives rise to the strong d"^ 
dependence in p^. It shoud be emphasized also that the Tr dependence is modified when 
the electronic states in the layers are not free Fermi gasses, and support longer-lived density 
fluctuations. As detailed in the following sections, this is the case, e.g., in the presence of 
disorder, and in the quantum Hall regime." 

IV. THE ROLE OF PHONONS 



The electron-electron scattering treatment predicts a transresistance po that vanishes 
at zero temperature as po r^ T^. This dependence is roughly satisfied in the experiments, 
confirming the dominance of the electron-electron interaction. However, as mentioned above, 
the experiment also shows a noticeable deviation of po/T'^ from a constant as a function of 
temperature, showing a maximum for T ~ 2K. For this samples the Fermi temperature is 
roughly QQK. The overall temperature dependence and the position of the maximum are 
very similar for different values of the layer separations rf, and the magnitude of pn/T"^ varies 
very little for three barrier thickness d = 500, 225 and 500A after the Coulomb scattering 
contribution was subtractedli^l. Also, for interlayer separations d = 500A, the observed pd is 
simply too big to be accounted for by the Coulomb interaction alone. This led Gramila et 
alx^ to propose an additional scattering mechanism. The obvious candidate: phonons. The 
phonon-mediated coupling between electrons in doped semiconductor layers separated by an 
insulating region of thickness ~ 100/im was studied by Huber and Shockleycil. Real phonons 
were found to be responsible for the interlayer interaction. The possibility of drag due to 
phonons was also proposed by Gurzhi and KopliovichB3. The first qualitative hint of the 
mechanism being phononic is the fact that the measured temperature dependence of Pd/T"^ 
resembles the acoustic-phonon-limited mobility scattering rate r^^ for two-dimensional 
electrons in GaAs. At high temperatures t~^ is linear in T, but crosses over to T^ or T^ in 

9 



the low-temperature Bloch Griiniesen regimetj, where the thermal phonon wave vectors are 
less than 2A;^. For the electron density of the samples in Ref.Ej the crossover occurs at a 
few degrees Kelvin. Thus, the temperature dependences of t^-^ /T"^ and po/T'^ are broadly 
similar. Furthermore, the dependence of p^/T^ on the relative electron density between the 
layers indicates that phonons could be playing a role. An electron in layer 1 decays through 
a backscattering process by emitting a phonon of wave-vector 2kF,i, and the phonon will 
transfer its momentum most efficiently to an electron in layer 2 if kp^i = kF,2- This implies 
that Pd/T^ should have a maximum when the electron densities are matched in the two 
layers, which is experimentally observedt^. Now, interactions of acoustic phonons with 
electrons are relatively weak in GaAs to account for the value of the observed transresistance. 
The proposed mechanism will be the exchange of virtual phononsEj'u, a process in which a 
phonon is emitted by one layer and then absorbed by the second without conserving energy 
from the electronic transitions. When the energy conservation constraint is relaxed, the 
phase space for scattering increases. Also, since the layers are separated by distances much 
smaller than the phonon mean free path, the phonons retain their phase coherence for the 
interaction between the layers. This two effects imply an enhancement in the transresistance 
due to virtual phonons. Tso et alB^ presented diagrammatic calculations including exchange 
of virtual phonons with a good agreement with the temperature dependence observed in the 
experiment (See Figure 2). 

The distinction between "real" and "virtual" phonons is not so clear cut if one extends the 
treatment leading to Eq. (^ to include a phonon-mediated Coulomb interaction. The force 
operator on a given layer will involve the phonon operators, and the force-force correlation 
function will include an electron-phonon interaction M{Q) and a the propagator of a phonon 
from layer one to layer two. The result is that one can write Eq. (P) with U^^\q, cu) replaced 
by an interaction T)i2{q, co) of the formcZl 



2^Q 



c<j2 — UJ% 



(17) 



^i2(g,^)=/|^|M(Q)|V«^'^ 

J Inh 
with ujq = VsQ the frequency of the acoustic phonon (omitting the distinction between 
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longitudinal and transverse) of wave vector Q = Jq^ + Q1. The term in square brackets in 
the above equation is the phonon Green's functions. The phonon mean free path Eph is not 
included. In the above expression, one can interpret huj as the energy transferred between 
layers and huoq as the energy of the intermediate phonon. When the denominator in the 
square bracket of Eq. (p^ vanishes, energy is conserved in the intermediate state. As 
pointed out by B0nsager et a/.EII, this expression contains both the real and virtual phonons: 
real phonons correspond to ImT'12, whereas virtual phonons correspond to Re 1^12. If we 
insert in Eq. ([T7| ) the long wave-length limited \M{Q)\^ = hQD'^/2pVs, we obtainE3 

9 

V,2{q, iv) ^ I '^,,, -- exp {-qd^l-uj^/[vsqy) , (18) 

q^l-uj^/{Vsqy ^ ^ 

indicating that the effective interaction diverges as g ^ ^jvs- Even though T>i2{q, uj) involves 
a small prefactor, if we substitute in Eq.(|^) the obtained expression for the coupling, the 
divergence \Vi2{q,u!)\'^ ~ |c<j — Vsq\^^ gives rise to a divergent pn- Although this is a 
spurious divergence that is removed by inclusion of a finite mean free path of the phonons or 
dynamical screening of the interaction, the large contribution explains partly why a "weak" 
phonon mediated interaction can compete with the Coulomb interaction as a mechanism 
for drag. For a phonon mean free path £ph below a critical value which for GaAs can be 
of the order of 0.2mm, the predicted distance dependence of the drag is ln(£ph/(i). The 
experimental distance dependence is yet to be clarified, although the evidence is to a weak 
dependence of the phonon contributionE^I. 

V. THE EFFECT OF DISORDER AND LOCALIZATION 

The effect of disorder was first studied by Zheng and MacDonaldt3, who included the 
fact that the density response function at small frequencies and small wave-vectors is given 
by 

, , dn Dq^ . 

dp Dq'^ — luj 

11 



with D being the diffusion constants given hy D = /^/2r, / being the mean free path, 
r = l/vp the scattering time, and dn/dfi the density of states . This formula for the 
density response function is vahd for q < l/l and uj < 1/r. Following the derivation of 
Pd ~ T^ presented in Appendix A, Zheng and MacDonald obtain that the low temperature 
dependence is modified to po ~ T^logT. This temperature dependence can be seen to 
result from the low frequency and low wave-vector behavior of Imxj(g, c<j) ~ ^/q'^ for 
q > [uj/Dy/'^. (Note that in the ballistic regime Imxj(g, cu) ~ ^/q-) The contributions 
to the integral in g for g < (co/DY'^ can therefore be neglected. We can obtain the low 
temperature behavior of po as followslBEHl: 

g2 /"I/t (^2 /"l/' dn g^ 

PD--n ^/^ / duj^r -/ — ~7i iT^T^logT. (20) 

In the above expression we have considered the low wave-vector contribution from the 
Coulomb potential as a constant U^^\q,uj) ~ e^/[kTp(kTFdY]. Note the change in distance 
dependence of pu from d""^ to c?"^. This is a consequence of having treated the Coulomb scat- 
tering as a screened interaction. The logarithmic term originates in the different spectrum 
of particle-hole excitations for a ballistic from a diffusive system and from the dimension- 
ality of the system. Both the spectrum of excitations and the dimensions (two in this case) 
conspire to give the logatrithmic term in the temperature dependence. At low temperatures 
one can ignore the contribution of Im x to the screened potential given by Equation (|^), and 
replace the interlayer interaction by 

with /cxF = 2'Ke'^dn/dp the single layer Thomas-Fermi screening wave vector. 

The case of strong disorder, where the localization length ^ the of states within each 
layer is of the order of the distance between layers was considered by ShimshoniE^I. For the 
case of Anderson insulators, the density response function within each layer is taken from 
the self consistent theory of localization of Wollhardt and Wolflea: 



dn Dq^ 

dp, Dq^ — [iuj + r(u;2 — ujI)] ' 



xiQ^ ^) = — r..2 L-^.,,/^.2 ^2^v (22) 
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with ujo a restoring frequency that incorporates the effects of locahzation. Within the scheme 
of WoUhardt and Wolfie, ujq is related to the locahzation length ^ through ^ = vf/{\/2uoq). 
The parameter ujq in the above expression is responsible for the vanishing quadratically of 
the conductivity a(u;) at low frequencies: 

Rea{uj) = e^\iMuj/q^)lmx{q,uj) ~ [D / {r'' uj^^)]uj'' . (23) 

Note that the logaritmic divergence at low q is now cut-off by 1/^ and the T^ dependence 
is retained. In the strongly localized regime, when the localization length is of the order of 
the inter-layer spacing, the distance dependence of po is modified. Screening is not effective 
in this regime, so we can evaluate p^ using the unscreened interaction f/q = lixe^e"'^'^ jq. In 
this regime, the low q dependence Imx(g, cu) will be given by 

dn ujc? dn 1 r,. , ^, 

Replacing the above expression in Equation (||) we obtairS, 

2 

327r e'^ynOn ^'"""^' \d 



Po^^AI^]it..dm. (25) 



If the localization length C, ^ d one has to include the effects of screening for wave- 
vectors q > l/{^'^kTF), and the distance dependence changes to ln(A;TF O/i^Tpd)^- When 
^ ^ oo we recover the diffusive result. Another case treated in Ref (e3) is that of the so 
called Efros-Shklovskii insulators, where Coulomb interactions are important and increase 
the conductivity at finite frequency in such a way that a{uj ^ 0) ~ tu as opposed to tu^. 
This has the consequence that the integration in uj is infrared divergent and is cut-off, at 
finite T by incoherent phonon processes. As a consequence the drag resistance diverges as 
zero temperature in the following way 



Pd oc f{d,^)i — \ exp 



%Y' 



with /(t/,0 = l/i^df for ^ < rf and f{d,^) = 1/^^ for ^ >(i. Also A;bTo = e^e^, with e 
the dielectric constant. 
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VI. PLASMON ENHANCEMENT 

The discussion of plasmon modes of Section |T| focuses on the small q dispersion at 
zero temperature. Since the drag resistance is given by an integral involving Imx in the 
integrand, at low temperature there is no contribution from the plasmons, since they appear 
in a region of the plane (g, cc') where Imx = 0. The function Imx "counts" the number of 
particle hole excitations of mementum a and frequency uo: it is nonzero in regions of the 
(g, uj) plane where particle-hole excitations are allowed. The plasmons are collective modes 
that are outside of the particle-hole continuum. For a more detailed discussion of this point 
the reader is referred to Nozieres and PinesE3. 

Flensberg and HuE3'E3 made the very interesting observation that at higher tempera- 
tures (T ~ Tp)-, Imx 7^ in the region of the plasmon pole, and therefore there will be a 
large contribution at intermediate values of q and one should observe an enhancement of 
the resistanceE2. Qualitatively one can understand this effect as a Schottky-like peak that 
develops from a thermally populated dissipation channel that is not available at zero tem- 
perature. In the plasmon-pole approximationl2Zl one approximates the dielectric constant 
e(g, cu) [which is the denominator in Eq. (H)] by 

e(g,^) = 2n(g)e-«'^|/?±(g)[o.-u;±(g)]+nmx(g,^±)|, (26) 

with I3± = [dRex{q,^)/d^]\uj=u}±- For small values of the imaginary part of x one can 
approximate the Lorenzian by 6 functions and obtain 

which we can replace in Eq. (^ to obtain the plasmon contributions to the drag rate: 

^ h^(3 /-"-i dq 3 Imx[q,uj±{q)] 

^"^ Txriin^e^AJo {2Tiy A\^^{q)\sinh\i3huj±{q)/2\ ^ ' 

The parameter q±^c defines the values of q for which the plasmon ceases to exist. The 
contribution to the drag from the above equation involves an integral of Ini x for frequencies 
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and wave vectors corresponding to the plasmon modes. From the above discussion, the 
contribution of the plasmons to p/j is zero at T = 0. At small temperatures, Imx at uj±{q) 
is small because the carriers do not have sufficient energy to be excited far above the Fermi 
surface. However, at intermediate temperatures (of the order of the Fermi temperature) 
there are enough thermally excited particles to give a large Imx at the plasmon poles and 
the drag will be dominated by this contribution. The numerical calculations indicate a 
maximum in po at a temperature of the order of 0.5Tp. For electron densities of the order 
of n = 1.5 X lO^^cm^^ for GaAs quantum wells, Tp ~ 50K. For T > 0.5 Tp the plasmon 
modes are no longer well defined since they can decay by emitting single particle excitations. 
This implies that the enhancement diminishes at high temperatures. 

The plasmon enhancement theory was put to test in the experiments by Hill et alxB with 
very good qualitative agreement. The temperature required to excite a plasmon appears to 
be lower than the value predicted by the theory, and the magnitude of the drag is larger 
than the prediction on most of the temperature range. This seems to imply that one should 
go beyond the Random Phase Approximation. The Hubbard approximation!^ provides a 
better fit to the data, but there is still room for improvement (see Figure 3). 

Improvements over RPA were considered by Swierkowski et a/.E2l in connection with drag 
experiments in electron-hole systems by Sivan et alM. The transresistance was calculated 
using Eq. (P) modifying the effective inter-layer interaction using the local field approach of 
Singwi, Tosi, Land and Sjolander (STLS)ElrE3. This approach includes correlations neglected 
in RPA, and the net effect is an increase of the effective interaction. This corrections use 
zero temperature local field corrections. Tso et alx3 studied a generalized RPA including 
exchange processes to infinite order in the Hartree-Fock potential and emphasized the fact 
that RPA is good only for very high densities. Thus, the RPA treatment of the Coulomb 
scattering overestimates the screening and renders the coupling weaker. The agreement 
with the experiments of Ref." is good at low temperatures. Giiven and TanatarE^ treated 
the plasmons and the phonons on equal footing using a the Hubbard approximation and 
found an enhancement due to the coupled plasmon-phonon modes. We still need better 
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approximations beyond RPA in the higher temperature regime to have a better quantitative 
understanding of the plasmon enhancement effects. 

VII. DRAG IN THE PRESENCE OF MAGNETIC FIELDS 

A. Hall drag 

In the presence of a magnetic field apphed perpendicular to the planes, there could in 
principle exist a trans-Hall resistance, that is, a voltage in the y direction of layer 2 when a 
current is applied in the x direction of layer 1. This drag Hall resistance is zero if computed 
it in the lowest order in the inter-layer interaction. The reason in that in lowest order the 
electron distribution in layer 2 is the equilibrium one, and also the fact that the momentum 
is being transferred from layer 1 in the direction of the current. Since there is no net current 
in layer 2 there is no Lorentz force and therefore no net Hall voltage. To see how this emerges 
formally so we can repeat the steps that lead to Equation p. In that case we obtain 



xy T^^P f d<i ,rr(s)f m2 f^ , Imxi(q,^)lmx2(q,^) 

7min2e^AJ (27r)^ " J-oo smh {(3huj/2) 



Pd 



where the functions xi(q, c<j) have to be evaluated in the presence of a magnetic field 
B(x) = V X A(xJ. The reason why the equation retains its form is the fact that the force 
operator is still given by Equation ( [A12| ) even in the presence of a magnetic field since 
[A(x,.),pq] = 0. In other words, the interlayer force does not depend on the applied field. 
Due to rotational invariance xi(q, tu) = xi(|q|,C(j) and the integral above vanishes by parity. 
The result p^ = was also shown by Kamenev and OregE3 using a diagrammatic approach. 
However, as emphasized by Hue3 this is not a general result. In particular, Hu showed that 
using an approach based on the Boltzman equation, if one includes an energy dependent 
lifetime r(e) one obtains a finite p^j^ . No experiment has reported so far a finite Hall drag, 
at least in the situation were one can guarantee that there is no tunneling between lay- 
ers. Patel et alx^ measured frictional drag in the presence of a field for modulation doped 
GaAs/AlGaAs double quantum wells separated by a barrier of lOOA for which tunneling 
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is significant. They found oscillatory behavior as a function of magnetic field in both the 
longitudinal and transverse drag resistivities. However, since for this interlayer separation 
tunneling is significantpS, our analysis does not apply to this case. 



B. Magneto Coulomb drag 

The physics of electrons in two dimensions under quantum conditions of temperature 
and magnetic field has been the subject of intense research in recent years, and the reader 
is referred to some of the very good reviews on the progress in understanding of the integer 
and fractional Quantum Hall Effectai^H. In this section we review the research on current 
drag in the presence of magnetic fields. 

In the presence of a magnetic field, the polarization function x(g, uj) assumes different 
forms depending on the number of filled Landau levels and on the strength of disorder. 
For example, for non-interacting electrons, if the density and magnetic fields are such that 
the outmost occupied Landau level is completely filled, the lowest energy for a particle- 
hole excitation is huc, with Uc = eB/mc the cyclotron frequency. Therefore Irachi will 
consist of a series of delta functions separated by Uc = eB/mc. If the outmost occupied 
Landau level is partially filled there will be excitations or zero frequency that correspond to 
transitions between orbits in the same Landau level. The presence of disorder smears this 
delta functions. One therefore expects that the transresistance will be sensitive to Landau 
level quantization for fields large enough that the cyclotron frequency is much larger than 
the disorder induced lifetime of the orbits. This problem was studied by B0nsager et alSJc3 
and Wuo and QinE3, who treated the individual layers as non-interacting electrons in the 
presence of disorder consisting of short range impurities. In the presence of disorder the 
one-particle density of states g{e) is broadened, and B0nsager et al. chose to substitute 
the comb of delta functions centered at e„ = (n + ^)huc for a Gaussian density of states as 
derived by GerhardtsE3 
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-2-^-^" 



(30) 



with £b = Jh/eB the magnetic length, and Fq = {2/'K)fiuJc{fi/T), r being the transport 
hfetime in the absence of magnetic field. The equations for the transresistance were solved 
numerically and revealed a very interesting "twin peak" structure as a function of filling 
fraction: as the filling factor is changed from an odd value (where the highest Landau level 
is half filled in the spin-unpolarized situation) towards the even value, the transresistance 
goes through a maximum before it is suppressed. The authors explain this non-monotonic 
effect as being the result of the competition between phase space for scattering and the 
strength of the effective interaction. When the Landau levels are completely filled, the 
chemical potential is in a gap, there is no dissipation, and the transresistance vanishes. This 
corresponds to the plateau regions of the Quantum Hall Effect. Also, the screening is 
strongly depressed in the region where the density of sates is small. In the transition region 
there are excitations of zero energy, and since density of states increases as one approaches 
the center of the band (half filled Landau levels) one gets an increase in x(g, u) contributing 
to an increase in p^. On the other hand, an increase in x(g, cu) means an enhancement of 
screening, or equivalently a decrease in the effective interaction. The prediction is therefore 
that p£, should roughly fulfill the relation 

Pd = gig2\Wi2\^, (31) 

with Qi the density of states at the Fermi energy for layer i, and W12 the effective inter-layer 
interaction. The striking twin-peak structure was observed in the experiments by Rubel et 
alxJ with good quantitative agreement with the theory. The twin-peak structure is observed 
up to a filling fraction u = 15. There is a marked increase in the value of the transresistance 
from p£, = SmQ at zero field and T = 3.1K to values of the order of 1 — 2^2, as expected 
from the theory (see Figure 4). The calculation does not include effect of localization and 
one should keep in mind that the Qi of Equation (^) refers to a density of extended states. 
The enhancement of po in the critical region was studied by Shimshoni and SondhiEa. They 



show that anomalously slow relaxation of density fluctuations at criticality yields a power- 
law p£) ~ 2"2-77^ where 77 is the anomalous diffusion exponent. 

The temperature dependence of po also shows different behavior with or without an 
applied magnetic field. In the discussion of Section |V| we found a low temperature de- 
pendence Pd ~ T^lnT in the diffusive regime. The range of temperatures for which this 
dependence applies corresponds to ksT ^ ^/r, with r a scattering time. This defines a dif- 
fusive temperature that for high mobility samples is of the order of 50 mK. Moreover, Zhen 
and MacDonaldt3 estimate-after including the different distance dependences in the ballistic 
regime and the diffusive regime-a crossover temperature of IQ^^'^^K ! At low temperatures 
for i? 7^ one expects the same T^ In T behavior since the motion is also diffusive, but 
the temperature scale where the effect sets in can be higher and experimentally accessible. 
The reason is that in the presence of a field one can think of diffusion as hopping between 
adjacent orbits of radius Cb *^ vpT. This implies that the diffusive form of the polarizability 
of Equation (|T9|) is now valid for q ^ 1/ ^b ■ The numerical solutions of B0onsager et al^ 
indicate that, for a choice of parameters corresponding to identical layers with a density 
n = 3 X 10^^771"^ and an interlayer distance d = 8OOA, the diffusive behavior sets in at 
T = 0.4K. 

Another very interesting aspect of the temperature dependence at finite magnetic fields 
predicted by B0onsager et a/.Ell and observed in the experiments by Rubel et a/.E3 is the 
dynamic screening originating a maximum in p^/T"^, the maximum being at Tp. The theo- 
retical reason is the same leading to plasmon enhancement discussed in Section |V^ . 

Wu et a/.Eilalso studied the interplay between Landau quantization and transresistivity by 
solving numerically Equation (P) in the presence of a magnetic field. They find pronounced 
oscillations as a function of field but not the "twin-peak" structure. Also, experiments by 
Hill et al!^ show oscillations in pD without the "twin-peak" structure. The oscillations are 
pronounced at low temperatures T ^ 4K. For high temperatures the Landau quantization 
effects are washed out and pd has a dependence 
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Pd oc TE^, 
which has not been addressed theoretically yet. 

C. The case of z^ = 1/2 

The case of Landau levels of filling fraction z/ = 1/2 (half filled) or even denominators is 
specials. The quantized Hall effect does not occur in this case, and the low energy physics 
is that of a Fermi liquid, behaving in many ways as electrons in zero magnetic fields. The 
theoretical approach that has proved most useful for understanding this state is the fermion 
Chern-Simons theoryEiJ, which is based in turn in the composite fermion theory developed 
by J. K. Jainil. 

For the fractional Hall effect, Orgad and LevitE3 and Duana studied the Coulomb drag for 
edge excitations using a Chern-Simons theory. The case of interlayer friction was considered 
for V = 1/2 by SakhiE3 , Ussishkin and SternH and Kim and Millig^. The dominant low 
temperature for p£, is found to be T^/^. This temperature dependence results from the slow 
diffusion of the density modes at filling fraction 1/2. In the composite fermion picture, at 
V = 1/2 the density response at small frequencies and small wave vectors is of the formcSl 

3 

with Xo = dn/dfi the electronic compressibility. The form of is similar to that of the diffusive 
regime at -B = of Equation (|TUp with an effective diffusion constant that vanishes linearly 
in q. This means that the long wavelength density fluctuations relax very slowly leading to 
an increase in the transresistance. To obtain the low temperature dependence we proceed in 
the same way as for the diffusive case, except that now, for small q, Imx ~ 1/q^, with the 
divergence now being cutoff at g ~ 2 {xo^kp) cu^/^ = kp.{uj/ujQy/^. The low temperature 
behavior is given by 



{kT^Ydf Jo e/3^^ + e-^^^ - 2 Jk^iu^/c^oV/^ q^ {k^Fd^ ' 
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The case of z/ = 1/2 poses some fundamental questions that are unresolved at the mo- 
ment of writing this review. The recent experiments by Lilly et alM for modulation doped 
GaAs/A^Gai_2:As double quantum wells separated by 200A, the drag resistance po has 
a qualitatively similar behavior as the longitudinal resistance of a single isolated layer p^x 
when the magnetic field is varied: where p^x is maximum due to the quantum Hall effect (for 
example at z/ = 1 and 2/3), so is po- As pointed out in Ref.E3 this is not surprising since both 
resistances are controlled by the density of states available for scattering. There is, however, 
a notable difference in their respective temperature dependences. Whereas pxx increases by 
only 6% as the temperature is lowered from T = 4K to 0.2i^,the drag resistance decreases 
by a factor of 40 in this temperature range. The numerical value of the drag resistance 
at -B = IIT, ( the field corresponding to z/ = 1/2) and at AK is about 2000 times larger 
than pd at zero field. In addition, the temperature dependence of pn differs qualitatively 
at u = 1/2 from the corresponding resistance at zero field: the experiment does not show 
evidence of a "phonon peak" in po/T"^, suggesting that phonons are relatively unimportant 
as a scattering mechanism in the z/ = 1/2 case. 

The most intriguing part of this story is the evidence of a finite value of pr, ~ 5^2 when 
extrapolated to zero temperature (see Figure 5). This is clearly in conflict with the above 
discussion leading to the temperature dependence of Eq. (|33|) that predicts a vanishing 
resistance at zero temperature. Since the scattering giving rise to p£, is inelastic, and the 
common view is that all inelastic processes cease to be effective at zero temperature, this 
experiment is inviting us to reanalyze the mechanisms of interlayer dissipation at zero tem- 
perature. This interesting experiment places current drag in the same arena as some recent 
efforts in understanding the issues of dephasingoS and resistance at zero temperatureEHl 
due to "inelastic" mechanisms. Some very recent papers address the u = 1/2 from different 
angles: Ussishkin and Sternly attribute the anomaly to pairing fluctuations whereas YangpJ 
suggests that the low temperature behaviour of the drag resistance is due to the interlayer 
distance being close to the critical value where the two layers form a collective incompresible 
state. 
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VIII. NON DISSIPATIVE DRAG 

The possibility of a drag effect at zero temperature was considered by Rojo and MahanO, 
who considered two coupled meosocopicEa rings that can individually sustain persistent 
currents. The mechanism giving rise to drag in a non-dissipative system is also based on 
the inter-ring or inter-layer Coulomb interaction, the difference with the dissipative case 
being the coupling between real or virtual interactions. One geometry in which this effect 
comes to life is two coUinear rings of perimeter L, with a Bohm-Aharonov $i flux threading 
only one of the rings (which we will call ring one). This is of course a difficult geometry 
to attain experimentally, but has the advantage of making the analysis more transparent. 
Two coplanar rings also show the same effectEHl. If the rings are uncoupled in the sense 
that the Coulomb interaction is zero between electrons in different rings, and the electrons 
are non-interacting within the rings, a persistent current Jq = —cdE/d^i = evp/L will 
circulate in ring oneo. If the Coulomb interaction between rings is turned on, the Coulomb 
interaction induces coherent charge fluctuations between the rings, and the net effect is that 
ring two acquires a finite persistent current. The magnitude of the persistent drag current 
Jd can be computed by treating the modification of the ground state energy in second order 

(2) 

perturbation theory AEq , and evaluating 

(34) 



dAE', 
-Id — - " 



d^2 

*2=0 

with $2 an auxiliary flux treading ring two that we remove after computing the above 
derivative. In other words, the persistent drag current is equivalent to a induced diamagnetic 
current that is finite even when the flux is zero in that system. (Note that a diamagnetic 
current is in general given by the differential change in energy with respect to a change in 
magnetic field.) The correction to the energy resembles the van-der Waals interaction, and 
its relevance to systems that can individually break time reversal symmetry was studied in 
Ref.Sn. The second order correction is given by 



^ep = - e i^.r r du r dj^-^^^t^^^^t!k^^ (35) 
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with Si{q,uj) the dynamical structure factor of ring i (see Eq. [A14| ). For a niesoscopic ring 



with an apphed Bohm-Ahronov flux one has to retain the discreteness of the spectrum: 

^.(^,a;) = ^r(g,..-|^g^), (36) 

with 5*1 (g, to) the structure factor at zero flux, (pi = $j/0O) aiid (pQ = hc/e the flux quantum. 
Note that, for a mesoscopic system, due to the presence of the flux, S^ (g, tu) ^ S^ {—q, uj). 
For the inter-ring interaction we take the unscreened Coulomb interaction (screening is not 
effective in one-dimension), so Uq = Ko{qd), with d the distance between rings and -^'0(2;) 
the modifled zeroth-order Bessel Function. In the limit of the interparticle distance much 
smaller than the distance between rings {kpd ^ 1), we obtain 

1 1 

^^ = '^\kpao)Hkpdr ^^^^ 

with oq. For dimensions corresponding to the experiments measuring persistent currentsEI 
J^) c^ 10~^Jq. The drag current is itself mesoscopic, and therefore vanishes in the limit 
of inflnite length. For a ring with a single channel carrying the current, an extension of 
the argument presented by VignaleEa for a bound in the value of the persistent currents 
gives Jd < Jq in general. The non-dissipative drag for two concentric rings was studied by 
Shahbazyan and Ulloa0 using a Luttinger liquid formulation of mesoscopic systemg53. They 
found that the inter-ring interaction modifles the period of the Aharonov-Bohm osillations. 
Related work on the effect of the interactions on the flux dependece was reported by Canali 
et alE. 

One basic difference between the dissipative drag in semiconductor systems and the non- 
dissipative drag in mesoscopic systems appears if one opens ring two so that no current can 
circulate, and computes and induced charge modulation, which will play the role of a drag 
voltage. The induced "voltage" is zero in this case. This can be seen by starting with a 
setup that in the absence of flux in system one is "parity even" . By this we mean that the 
charge distribution in wire two is symmetric around the center. Let us call T and V the time 
reversal and parity operators respectively. We want to compute the induced dipole moment 
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0, (38) 



in ring two X2 = (\l/o|5;2|^I/o)- The coordinate operator satisfies VTx2{VT)~^ = — X2 while 
the wave function is invariant under VT, which imphes x = 0. The charge distribution in 
ring one remains uniform, and there is no induced vohage. 

A natural candidate for the study of non-dissipative drag other than mesoscopic systems 
are superconductors, which can sustain macroscopic persistent currents. The extension to 
this case was done by Duan and YipE3. For wires the corrections to the zero point energy 
of the charge fluctuations (plasmon modes) was computed when two superconducting wires 
individually carry a supercurrent. In analogy with the discussion of Section || the dispersion 
of the coupled modes is given by the determinantal equation 

{u — qviY — s'^q^ —X 

—X {uj — qviY — s^q^ 

with s the plasmon velocity of a single wire, X = A7mQe'^q'^KQ{qd)/m comes from the inter- 
wire interaction, and f 1^2 are the superconducting velocities of each wire. For the case of 
wires the two couple modes of frequencies uj±{q) are linear in q. The zero point energy Eq = 
J2q^ [^^+(?) + ^^~{q)] depends now on the relative superfiuid velocities. The superfiuid 
velocity is /i(V\l' — 2eA/c)/2m with \1' the order parameter. The supercurrent in wire two 
I2 is computed in an analogous way as Eq. (^4|) by taking the derivative of the free energy 
with respect to the vector potential in wire 2, and 

h = —{P22V2 + P21V1), 
m 

with P12 = hri^e^ / IQTmi^ s^ d? representing the drag term. If one starts with a situation in 
which there is no current in either wire and slowly increases the current in wire one, the 
prediction is that a current will start to flow in wire two of magnitude ep2iVi/m. Supercon- 
ductivity is essential to this effect, so that the wire can be trapped in a metastable state. 
Extensions of the above arguments to superfiuid Bose systems were presented by Shevchenko 
and Terent'evS and Tanatar and DasE^. The case of the transresistance of an excitonic con- 
densate with electrons in one layer and holes in the other layer was studied by Vignale and 
MacDonald^a, who found a discontinuous jump in p^ at the condensation temperature. 
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Two groups attempted to measure the non-dissipative drag. Giordano and Monnieii^ 
measured the drag between a superconductor (Al) and a normal metal (Sb), and found a 
non-reciprocal drag effect that was finite only in the transition region, that is, at temper- 
atures close to the critical temperature of the Al layer. Similar results were reported by 
Huang et alM with the superconducting system being AlOj; and Au for the normal metal. 
An interpretation of the non-reciprocity effect in terms of inductive coupling of the spon- 
taneously generated vortices in the superconductor and the normal metal was proposed in 
Ref.E3. It seems clear that these experiments do not provide evidence of the supercurrent 
drag. 



IX. SUMMARY 

We have reviewed the recent theoretical and experimental efforts in the growing field 
of current drag. For electronic systems separated far enough so that there is no tunneling 
between them, we found that through studies of the transresistance po one can extract 
information not only about the direct Coulomb interaction between electrons in different 
systems, but also about the collective modes of the coupled systems, and about phonons 
(virtual and real) that can propagate through the barrier separating the systems in question. 
The magnetic field B has a non-trivial effect on the drag resistance. The value oi p^ is larger 
af finite fields than its value at zero field, and shows a large enhancement in the transition 
region between quantum Hall plateaus. This effect is understood qualitatively in terms of 
an interplay between changes in phase space available for scattering and variations in the 
effective interaction with B. Some theoretical considerations, like those in connection with 
the non-dissipative drag, and the ones for strongly disordered systems, are still waiting 
experimental verification. The recent experiment on the z/ = 1/2 case pose a fundamental 
question connected to the general theory of transport: is it possible to have an intrinsic 
resistance at zero temperature due to the electron-electron interaction? The magnitude of 
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the challenge seems as big as that of the progress made so far in the field. 

We acknowledge interactions with W. Dietsche, J. P. Eisenstein, A. J. Leggett, A. H. 
MacDonald, G. D. Mahan, C Proetto, E. Shimshoni, S. UUoa, P. Vasilopoulos and G. 
Vignale. 

APPENDIX A: 

In this section we follow the derivation of Zheng and MacDonalclll3 of the interplanar 
resistance using the memory function formalism. The formalism of MoriS uses the projector 
technique to write the response function in terms of a memory function. A simple case is 
a classical particle subject to stochastic forces, for which one is interested in the velocity- 
velocity correlation function (f){t) = {v{t)v{0)). The Laplace transform is 



with K{s) the Laplace transform of the memory function 



(Al) 



Here F{t) is the force acting on the particle, and m(f^(0)) = /c^T. The friction coefficient 
is then given by a force-force correlation function. 

For the quantum mechanical case, the calculation of an inter-layer resistance is then 
performed by identifying the relevant memory function. The starting point is Kubo's formula 
for the conductivity: 



a,,{uj) = ^J^ dte'-'{J,\Mt)) (A3) 

where (3 = l/ksT, A is the area of each of the two-dimensional systems considered, i and j 
are the layer indices, and J is the zero wave vector component of the total current operator. 
The inner product appearing in (|A3|) is 
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{Jj\Ji{t)) = \ I dXTi Uexp(Aif)J,exp(-AiJ)J,(t)l , (A4) 

p Jo ^ J 

with Ji{t) = eE'^Vje^E'^* = e^^^Ji, and C is the Liouville superoperator defined through 
its action on an operator O as CO = ^[H,0]. The next step is to use Mori's projector 
method to write the inverse matrix a^J(uj) in terms of a memory function, and hence obtain 
transresistance. Zhen and MacDonald define a superoperator V that projects onto the 
current: 

with {Ji\Jj) = Sij^Xi, and Xi = riie^/m. 

With this definitions we have (for Imuj > 0), 



a,,{u;) = ^ / dte'-'{J,\e-'^V,) = ^{J^\ 7\Jj) (A6) 



'ijv^; — "jy rfte*"^ (Jj|e * \Jj) — ~j{M _ ^I'^j/ 

= ^{J^\ ^! ^J J,) (A7) 

A uj — QC A uj — C uj — Ql 

' 00 f^^ u Axk ' 00- QC ' ^^' ^ ' 

where we have used VJi = iVC Ji = 0, which results, in the present case, from the fact 
that [Ji, Jj] = 0, meaning that the derived expressions are also valid in the presence of a 
magnetic field that breaks time reversal invariance. We are interested in the real part of 
a^^, which gives the resistance. From the above equations we have 

Pi2{oo) = ^^-Re{j,\^-\j,) (AlO) 

^XiX2 uj-QC 



(3 



1 fOO 

/ dte'^'{F^\e-'^^'\F2), (All) 

noc'^ Jo 



Anin2i 

where the force operator is Fi = —m/eJi. The contribution to the force due to the inter-layer 
interaction potential f/q is 

Fi = -F2 = ^^Eqf^.piVt (A12) 
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The first approximation is to replace e~^^^ by e~*^. The leading order in f/q corresponds 
to the correlation function in ( |A11| ) evaluated in the uncoupled case. Using a representation 
in terms of exact eigenstates, one finds (for the static limit) 

■nnin2e^A J [2'kY J-oc smh {(3huj/2) 

with Xj(g, tu) the density-density response functionQ of layer i: 

Imx.(g,c.) = ^(1 - e^-^)l Ee-'3^™Kn|p,(g)|m)|2<5[c. - {E^ - E^)/h] (A14) 

= (l-e^-^)5(g,^), (A15) 

and Zi the partition function of the isolated layer. 

Here we show that pu vanishes as T^. For simplicity we treat the layers as identical. 
The approximation used in the literatureEj'EO is to replace Xj(g, cj) by its non-interacting, 
zero temperature valuetO: 



Imx(g,^) 



Q{kF — \x_\)Jkp — x^_ — Qikp — \xj^\)Jk'p — x 



(A16) 



with x± = rauj/hq ± q/2. At low temperatures, the factor l/sinh^(/3^u;/2) in (H) ensures 
that the values of uj are small. Then we can separate the density response functions as 

, , 2m^id(2kF - q) ^, . /,,^,^ 

Imx g, uj) = , / =c^ = F{q)uj. A17 

nh^qJ{2kp)^ - g2 



^ / ^^ ■ u2/^fe /ON = (^^^) ^I3^(?) • (A18) 



The integral over uj is 

-oo~^^sinh2(/?;icj/2) ''"""^ ' 3/i 
The remaining integral is 

I = Jdqq'U^iq)F{qf. (A19) 

Here, one can replace the bare interlayer interaction with 

97rp2 p-Qd 

UM) - U^'^^ = ^^^, (A20) 

q eooe(g) 
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where e{q) is the effective dielectric function for the two parallel conducting planes, and d 
is the distance between planes. In the random-phase approximation, 

e{q) = 1 - 2U,{q)x{q, 0) + U^{q)x\q, 0)(1 - e-"^'). (A21) 

In the long wavelength limit Uii{q)x{q, 0) = — Q'tf/q'j with gxF = 2/ao the screening wave 
vector in two dimensions. The last term in e(g) dominates whenever gxF ^ 1- In this limit, 
the integral / can be approximated as 

9'^?7^-z^;;;7vJ (A22) 



3m4C(3) 



2el,h'klq^^d^- 
Collecting this integrals we obtain 



(A23) 
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FIGURES 
FIG. 1. Schematic representation of the memory function expression for the transresistance 

to lowest order in the interlayer interaction [Eq.(g)] [reproduced from Ref. (113)] 

FIG. 2. The scattering rate due to the Coulomb scattering and virtual phonons t^ /T^ as 
a function of temperature for different separations. Note that pn cc t^ . The solid circles are 
the experimental results of Ref. (□) and the solid lines are the theoretical results from Ref. (E3). 
Inset: Contribution to pD oc r^ due to exchange of virtual phonons as a function of temperature. 
[Reproduced from Ref. (E3)] . 

FIG. 3. The scaled transresistivity pt/T~'^ {pt = pn) versus the reduced temperature for 
different densities (the densities in both layers are the same) . The dashed (solid) lines are the RPA 
(Hubbard) calculations of Flensberg and Hutll, and the circles are the experimental results of Ref. 
(e3). [Reproduced from Ref. (Ej).] 

FIG. 4. The transresistance Rt (pd) ^^ ^ function of magnetic field B for a coupled electron 
gas with a separation barrier of 30 nm. shown for different temperatures (plotted with offset for 
clarity). The electron density in both layers is n = 3.2 x lO^^cm^^. The longitudinal resistance is 
also shown. [Reproduced from Ref. (^).] 

FIG. 5. Measured temperature dependence of pn at u = 1/2 (solid line). The broken lines 
are calculations from Ref.E3'cl of po assuming two different values of the composite fermions mass 
(dotted, m* = 12mf,, dashed, m* = 4m.6, where mf, is the GaAs band mass). [Reproduced from 
Ref. (S).] 
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